Semi-Analytical Calculation of the Rouse Dynamics of Randomly Branched Polymers 
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We present a semi-analytical approach to the determination of the dynamic properties of randomly 
branched polymers under the Rouse approximation. The principle procedure is based on examining 
a spectrum of eigenvalues which represents the average dynamic behavior of various structures. The 
calculated spectra show that the eigenvalue distribution is random even within a single structure 
which in turn produces a continuous spectrum of values for the entire class. The auto-correlation 
functions for the radius of gyration squared were calculated based on these spectra, which confirms 
that the dynamics is non-exponential as earlier reported. A universal stretched exponent is also 
found in this study. 
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Among macromolecular structures found in many im- 
portant systems such as plastics, proteins, and sol-gel 
networks ^J-|!J, randomly branched structures occupy a 
unique place, as they display distinct physical proper- 
ties from those of linear or regularly branched polymers. 
A typical molecule is constructed from a connection of 
linear polymer portions of various lengths at branching 
sites selected randomly according to the underlying phys- 
ical branching mechanism [see Fig. |. The study of 
randomly branched polymers (RBPs) poses an interest- 
ing challenge to condensed matter theorists, as the aver- 
age over many different branching structures must be in- 
cluded in addition to the ensemble average taken over the 
configurational space for a given single structure. Most of 
the previous studies on these polymers have been focused 
on their conformational properties , which reveal dif- 
ferent scaling properties from linear polymer chains. The 
more intriguing aspect of randomly branched polymers 
is probably their connection to the spin glass problem: 
there exists a difference in quenched and annealed dis- 
orders, typical to random systems The dynamic 
properties of these polymers, however, remains largely 
unexplored. 

In a recent study 0, we were able to show, with the 
aid of Monte Carlo computer simulations, that the dy- 
namic behavior of these polymers in a good solvent does 
not exhibit the conventional exponential relaxation nor- 
mally seen in other regularly branched structures. The 
difficulty associated with the dynamic study of this class 
of polymer is that each polymer has a unique dynamic 
behavior dependent on its structure. Although the struc- 
tural averages can be performed using various techniques 
developed in statistical physics Q, it is unclear how to 
construct a rigorous analytical theory to explain the dy- 
namic properties even for the simplest cases, which have 
been examined by using scaling arguments |lj|7j] . In prin- 
ciple, numerical simulations of polymers provide an al- 
ternative approach for examining the basic dynamic be- 
havior of these randomly branched structures, allowing 
the sampling of many different branching structures. In 




FIG. 1. Sketch of a typical randomly branched polymer. 

practice, however, the number of branching structures 
that can be sampled is limited, due to the extensive com- 
putational time required. 

There are many unaddressed questions regarding the 
dynamics of RBPs. For example, one simple question is 
whether or not the non-exponential behavior is indeed 
an intrinsic property of the randomness. If it is, then 
it will show up inherently even in the simplest type of 
polymer dynamics, the Rouse dynamics, where the ex- 
cluded volume and hydrodynamic interactions are both 
neglected. In this article, we examine the Rouse dynam- 
ics of an autocorrelation function and intrinsic viscosity 
of RBPs using a semi-analytical method with a satisfac- 
tory average over the structural space. 

For a given structure, the Rouse dynamics are charac- 
terized by a NxN Rouse force matrix A appearing in the 
Langevin equation |10|, 



C^ = kA?+f(t) 



(1) 



where £ is the inverse mobility, r is a N-dimensional vec- 
tor containing the positional coordinates of the N 
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FIG. 2. The average eigenvalue spectrum for RBPs of 
N = 100. The spectrum appears continuous suggesting that 
the eigenvalues are distributed randomly for a single struc- 
ture. The inset graph shows dominating eigenvalues in the 
long time scale, and emphasizes the continuous distribution 
of the modes. The five peaks which extend beyond the range 
of the main graph, reach values of approximately 0.7, 2.5, 0.7, 
0.1, and 0.03, from left to right. 

monomers, and f(t) is a iV-dimensional random force 
vector containing Brownian forces acting on these N 
monomers. The N eigenvalues of A determine the relax- 
ation rate and dynamic properties of various collective 
normal modes of the given structure. 

The matrix A is directly determined by how the struc- 
ture (Fig. 0) is labeled and connected, which is ana- 
lytically generated here by a cut and paste algorithm, 
similar to the one used by Cui and Chen H. Starting 
from a linear polymer, we choose a bond randomly and 
cut the molecule into two portions. The smaller portion 
was then attached to a randomly chosen monomer on the 
larger portion. The move was discarded if the anticipated 
bonding site was already a branching monomer. The spa- 
tial positions of the monomers were not kept track of as 
the Rouse matrix is only concerned with the connectivity 
of the structure when the monomers have no volume. The 
process was initially repeated 10 5 times to create an equi- 
librate structure before data collection. This algorithm 
was originally suggested by Rensberg and Madras ||, in 
their study of lattice "animals" , and has been shown to 
yield a branching ratio between the number of branching 
nodes, n§, and total number of bonds, N, of n^/N r* 0.25 

A Rouse matrix was created after every ten cut-paste 
moves, to ensure that the structures were not all corre- 
lated. Once the Rouse matrix was defined, a complete 
set of eigenvalues were then determined numerically 0, 
and binned to create a normalized histogram. Each bin 
was given a width of 1.0 x 10~ 4 to create a fine scale im- 
age of the distribution. For each given N, a total of 10 7 
structures were examined. 

The spectrum that is produced is distinctly different 
from that of a linear polymer. The spectrum for a linear 



polymer would be constructed of only N eigenvalues, and 
thus would appear as N discrete lines spaced according 
to the following equation pi , 



A„ = Ai p 



P 



1,2, 



(2) 



where Ai is the smallest eigenvalue. In Fig. g the eigen- 
value histogram for a RBP structure of N = 100 bonds 
is displayed. The magnified portion for eigenvalues be- 
tween and 1, which governs the longest time scales of 
relaxation, shows that the distribution is continuous, and 
roughly follows a power law near A = 0. This histogram 
distribution indicates that the eigenvalues are uniformly 
distributed due to various segmental lengths and the cou- 
pling of the segments in the molecules. The figure also 
demonstrates that the eigenvalues of a single structure 
are randomly distributed, opposed to an ordered distri- 
bution as in a linear polymer. These features can be seen 
more clearly in Fig. |, where the histogram for N=10 is 
displayed. In a structure this small there are a limited 
number of distinct structures, and thus a limited number 
of possible eigenvalues, so the eigenvalues appear discrete 
in the spectrum. The figure clearly shows how the eigen- 
values are randomly arranged in the spectrum. 

Returning to Fig. ||, there appear to be several eigen- 
values which occur more frequently than the average. 
These eigenvalues are most likely caused by frequently 
occurring sub-structures that are regularly branched or 
linear. For example, the relative motion of a sub-unit 
that has two outer monomers connected to a stationary 
branching unit, will have an eigenvalue of 1 |Tl] ]; the large 
peak indicates the frequent appearance of such units. 

One of the most direct probes of the internal structure 
is the autocorrelation function associated with the radius 
of gyration squared. As has been appreciated in other 




FIG. 3. The average eigenvalue spectrum for RBPs of 
TV = 10. The small size of the polymer generates a distinctly 
different spectrum from that in Fig. ^| The spectrum is dis- 
crete, yet randomly distributed. This spectrum still generates 
a non-exponential behavior. 
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FIG. 4. Auto-correlation functions for a linear polymer (a) 
and a RBP (b) of length N=100. The function in (a) has 
been scaled such that the largest characteristic time, fi, is 
set to match that in figure (b). The other largest character- 
istic times are also placed on the graph to demonstrate their 
magnitudes. In figure (b), long time relaxation is dominated 
by multiple modes, which cause the observed non-exponential 
behavior. 

macromolecular structures, such a correlation function 
characterizes the radial relaxation motion, and defines 
the so-called elastic relaxation for a macromolecule [ [f2[ . 
The dynamic relaxation of such a correlation is comprised 
of the various internal relaxation times. Due to the large 
number of small eigenvalues in a RBP, no eigenvalue can 
be singled out to yield dominating exponential behavior. 
More specifically, we examine 
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(3) 



where t is the rescaled time such that i — 2kt/(, and 
iJ(Ai), represents the probability of the occurrence of 
the eigenvalue Xi given by the histogram. 

To further illustrate the significance of many time 
scales involved in the dynamics, the autocorrelation func- 
tion for the radius of gyration squared (Eq.||) was calcu- 
lated for a linear polymer and a typical RBP of N=100 
monomers in Fig. [|a and Fig. ||b respectively. Also 
plotted in the figures are the longest characteristic time 
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FIG. 5. A semi-log plot of the correlation function vs. 
rescaled time, for the RBP in figure [|b. The curve is fit- 
ted with an exponential function to demonstrate the sys- 
tematic deviation of the curve from this function, implying 
a non-exponential relaxational behavior. 

scales, indicated by the vertical arrows. The RBP has 
a significantly shorter relaxation time than the linear 
polymer since the structure is more compact, thus in 
the comparison of the distribution of the characteristic 
times, the linear polymer's time scales have been scaled 
such that the longest characteristic times of both poly- 
mers match exactly. Of course, when t is extremely large 
only the largest time scale is dominant, but this only 
occurs when the correlation is no longer significant, thus 
the correlation function is dominated by non-exponential 
relaxation. Comparing the two figures, we see that the 
characteristic time scales in a linear polymer, in which 
t p = Tip~ 2 ,p = 1,2, ... ,N [[To), are sparsely spaced in 
the region of interest, while the time scales in the RBP 
are irregularly distributed in the same regime. The two or 
three longest time scales in Fig. [I|b, are sufficiently close 
together to compete with each other; the net result is a 
multi-exponential function rather than a single exponen- 
tial function as in Fig. ||a. When averaged over different 
structures, various long time scales must be taken 
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FIG. 6. A double logarithmic plot of the correlation func- 
tions, for N=10,25,50,75,and 100. These plots are used to 
calculate the parameters a and f in Eq. (W). 
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into account, which yields an average over various multi- 
exponential functions. As discussed by Palmer et al. [fl3|| , 
and Iori et al. pH , such a process produces much slower 
decaying correlation functions, which sometimes are rep- 
resented by various non-exponential functions. 

The co-operative behavior of the competing time scales 
can be seen in Fig. ^, from another perspective. The fig- 
ure shows a semi-log plot of autocorrelation function for 
the same function in Fig. ||b. The dashed line is the best 
fit to the solid curve with an exponential function, which 
demonstrates that the correlation function is not expo- 
nential as in linear polymers. This result is in agreement 
with previous Monte Carlo simulation results which show 
the same non-exponential relaxation Qj. 

We have shown thus far that the eigenvalues of RBPs 
are randomly distributed, and the continuous distribu- 
tion near A = yields insight into the non-exponential 
relaxation. Without further knowledge of an analytical 
expression of H(X), the summation in Eq.(||) cannot be 
carried out exactly. As in previous work |Q, due to sim- 
ilarities of this system with spin glass systems and other 
random disordered systems |l3|-|l6||, it is instructive to 
compare the correlation curve to a stretched exponen- 
tial. 

C st (i) = (4) 

where f represents a general characteristic time for the 
entire class of RBPs, and a is the stretched exponent. 
As can be inferred from Fig. [| the stretched exponen- 
tial function gives a rather good approximation with the 
curves produced from Eq. (0). The a values for differ- 
ent N are shown in Table K which demonstrates that 
a m 0.87 ± 0.01 appears universal for all large N consid- 
ered here. The range used to calculate the exponents was 
selected by choosing the range of time scales such that 
the auto-correlation function was approximately 0.8 to 
less than 10 -4 , corresponding to the region where an ex- 
ponent of one can be produced for linear polymers. This 
value can be compared with a « 0.8 exponent calculated 
in earlier work |tJ ; the discrepancy is not surprising since 
this model considers no excluded volume effects. 

To further determine the universality of the stretching 
exponent a, we have calculated the eigenvalue spectra for 
RBPs with a different branching probability, correspond- 
ing to 113/ N = 1/2. In these polymers there exist only 
end and branching points in the structure, and no linear 

TABLE I. Stretched exponents, relaxation times and in- 
trinsic viscosity determined for various N. These systems 
correspond to a branching ratio of n^/N ~ 0.25. 



iV 


a 


T 


m 


10 


0.902 ± 0.005 


7.55 ±0.05 


1.49 ±0.05 


25 


0.879 ±0.005 


30.2 ±0.05 


2.81 ±0.05 


50 


0.871 ± 0.005 


87.8 ±0.05 


4.38 ±0.05 


75 


0.865 ± 0.005 


163.6 ±0.05 


5.63 ±0.05 


100 


0.861 ±0.005 


254.2 ±0.05 


6.69 ±0.05 
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FIG. 7. A log-log plot of f vs. N. Circles and squares 
correspond to two different classes of RBPs, those with 
riz/N = 0.25 and n^/N — 0.5 respectively. 

parts are allowed. To create these structures, a base star 
nucleus of 4 monomers was used. Using this base an 
end monomer was randomly selected and two monomers 
were attached to change the end monomer to a branch- 
ing point. This process was repeated until the desired 
structure size, N — 9, 25, 49, 75, and 99, was reached. 
The resulting correlation curves demonstrate the same 
stretched exponent as in the case of much lower branch- 
ing probability. The results are shown in Table |n]. This 
suggests that there is little or no dependence of the 
stretching exponent on the branching probability. Also, 
the examination of the exponents suggests that there 
is an insignificant correlation with N, suggesting that 
this stretching behavior is universal to the entire class of 
RBPs. Whether the quenched and annealed universal- 
ity classes |||| will produce the same exponents cannot 
be determined without the introduction of the excluded 
volume. 

The exponent f , displayed in Tables | and ||, in equa- 
tion ^ can also be calculated from Fig. ^. The plot 
against N in FigJ^ yields the scaling laws, 

f ~ at(!-5±o.i) j- Qr n ^ N _ 25 ( 5 ) 

t ~ 7v (1 - 2±al) for n 3 /N « 0.5 (6) 

Our early work suggests a larger exponent with the 
inclusion of an excluded volume interaction. This is con- 
sistent with the fact that the excluded volume effect slows 
down the dynamics generally |L7|. 

TABLE II. Stretched exponents, relaxation times, and in- 
trinsic viscosity for various N. These systems correspond to 
a branching ratio of nz/N ~ 0.5. 



N 


a 


f 


M 


9 


0.924 ±0.005 


4.4 ±0.05 


1.20 ±0.05 


25 


0.874 ±0.005 


17.3 ±0.05 


2.24 ±0.05 


49 


0.859 ±0.005 


39.0 ±0.05 


3.15 ±0.05 


75 


0.851 ±0.005 


63.0 ±0.05 


3.80 ±0.05 


99 


0.846 ±0.005 


85.1 ±0.05 


4.25 ±0.05 
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FIG. 8. A log-log plot of the intrinsic viscosity vs. TV. Cir- 
cles and squares correspond to two different classes of RBPs, 
those with n^/N = 0.25 and nz/N = 0.5 respectively. The 
straight line represents Eq.(Q). 

The intrinsic viscosity, [rj], is also a significant probe 
of the internal dynamics of a macromolecule, and can be 
directly measured experimentally One can show that [rj] 
is related to the eigenvalues |18| by 



lv} = 



1 



N A a 2 ( (N + 1) 



N 

E 



g(Ai) 
A; 



(7) 



where Na is Avogadro's number, rj s is the viscosity of 
the solvent, p is the weight per monomer, and a is the 
intermolecular bonding distance. The rescaled intrinsic 
viscosities are displayed in Tables | and |l[ The log-log 
plot of [rj] vs. N shows 



iv 0.63±0.01 for ns/N 
jy0.47i0.01 for n;}/N 



0.25 
0.5 



(8) 
(9) 



The exponent calculated here is for a Rouse model dy- 
namics. It is known that under the Rouse approximation, 
no simple scaling relation between [rj] and N can be de- 
duced. While linear polymers show rj oc N , star-burst 
dendrimers show r\ tx InN pd| . 

In practical good solvents, both hydrodynamic and 
excluded-volume effects must be considered. In such 
cases, a simple relationship can be deduced based on the 
fact that [77] is inversely proportional to the particle den- 
sity S 3 /N, where S is the radius of gyration. Based on 
this scaling argument, Daoud et al. [jlj have concluded 
a scaling exponent of 3/8 for [77], which, of course, can- 
not be directly compared with the exponent found here 
due to different physical origins involved. The expo- 
nent 3/8 agrees with the experimental results of branched 
polyethylene by Patton et al. ||. A more rigorous treat- 
ment of [77] for RBPs beyond the scaling argument, how- 
ever, cannot be found in previous literatures. 

In summary, by examining the average eigenvalue spec- 
trum of RBPs, we were able to show that the eigenvalues 
are randomly distributed as opposed to periodically orga- 
nized. The result of this distribution is the competition of 



relaxational modes which leads to a non-exponential be- 
havior in the relaxation of the polymer. Fitting the auto- 
correlation functions of the radius of gyration squared to 
a stretched exponential we found a stretching exponent 
of a = 0.87 ±0.04. 
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